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Harmonic Content of Strain-induced Potential Modulation in Unidirectional Lateral 

Superlattices 
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Detailed analysis of the commensurability oscillation (CO) has been performed on unidirec- 
tional lateral superlattices with periods ranging from a— 92 to 184 nm. Fourier analysis reveals 
the second (and the third) harmonics along with the fundamental oscillation for a>138 nm 
(184 nm) at low-enough temperature, evincing the presence of corresponding harmonics in the 
profile of the potential modulation. The harmonics manifest themselves in CO with demagnified 
amplitude due to the low-pass filtering action of the thermal damping factor; with a suitable 
consideration of the damping effect, the harmonics of the modulation potential are found to 
have the amplitudes V2 and V3 up to roughly 30% of that of the fundamental component Vi, 
despite the small ratio of the period a to the depth d~99 nm of the two-dimensional electron 
gas (2DEG) from the surface. The dependence of Ki on a indicates that the fundamental com- 
ponent originates at the surface, while the higher harmonics arise from the effect of the strain 
that penetrates down into subsurface. The manipulation of high harmonics thus provides a 
useful technique to introduce small length-scale modulation into high-mobility 2DEGs located 
deep inside the wafer. 

KEYWORDS: lateral superlattice, commensurability oscillation, harmonics, strain, stress, piezoelectric ef- 
fect, two-dimensional electron gas, GaAs, AIGaAs 



1. Introduction 

A high-mobility two-dimensional electron gas (2DEG) 
formed in a GaAs/AlGaAs-based heterostructure, often 
quoted also as modulation-doped field-effect transistor 
(MODFET), represents arguably the cleanest electron 
systems in solid state materials that have ever been 
made. The record mobility has reached a value as high as 
/i=3,100 m 2 /Vs, 1 which amounts to the electron mean- 
free path exceeding 0.1 mm. A MODFET, combined with 
modern nano- fabrication technologies, forms the basis 
of a wide range of experimental studies on the behav- 
ior of electrons under a designed environment. 2 In the 
majority of cases, the electrons in 2DEG are controlled 
by a patterned gate placed on the surface of the 2DEG 
wafer. Sizable efforts have been directed toward making 
the length scale thus introduced smaller, in pursuit of 
phenomena observable only when the artificial structure 
has a small-enough length scale. 3 A serious drawback of 
MODFET in this respect is the fact that the 2DEG re- 
side at the depth of tens or hundreds of nanometers from 
the surface; this is inevitable since the crux of attaining 
high mobility in a MODFET is to set back sources of 
electron-scattering, the ionized donor dopant (Si) being 
the most operative, as far away from the 2DEG channel 
as possible. The large distance of 2DEG from the surface 
renders it a formidable task to introduce a length scale 
smaller than the depth using nano-fabricated gates on 
the surface, since effects originating at the surface are 
generally expected to decay exponentially with the ra- 
tio of the depth d to the lateral length scale ai atora i as 

(X exp(-27rd/ai a teral)- 

A unidirectional lateral superlattice (ULSL) typifies 
the combined system of nano-fabricated gates and MOD- 



FET, where one-dimensionally modulated potential is in- 
troduced to 2DEG by a grating placed on the surface. 
ULSLs offer a unique opportunity to evaluate, through 
the amplitudes of commensurability oscillation (CO) ap- 
pearing in the low-field magnetoresistance, 4 the magni- 
tude of the potential modulation seen by the electrons. 5-7 
In the present paper, we investigate, by examining the 
CO in ULSLs, the harmonic content of the modulated po- 
tential introduced by elastic strain arising from differen- 
tial contraction between the gate material and GaAs. 8-14 
The period a of ULSL samples are chosen to be relatively 
small so as to be close to the depth d~99 nm. Although 
the fundamental component resulting from the modula- 
tion period a equal to that of the grating dominates the 
CO, second and third harmonics corresponding to pe- 
riods a/2 or a/3 are also detected, more clearly when 
the temperature is lowered. We discuss the importance 
of the thermal damping on the detectability in CO of 
the smaller period modulation. Modulation amplitudes 
V n (n=l, 2, 3) for the fundamental component and har- 
monics are explored as a function of a, which hints at 
the origin of each component. The harmonics turn out to 
be more persistent than simple exponential decay men- 
tioned above. Noticeably, a/2 and a/3 are significantly 
smaller than d. Therefore exploitation of high harmon- 
ics proves to be a powerful tool for introducing length 
scale smaller than the depth into the high-quality deeply 
residing 2DEGs. 

2. Experimental 

Five ULSL samples with differing periods (a=92, 
115, 138, 161, and 184 nm) were prepared from the 
same Alo.33Gao.e7As/GaAs single-heterostructure 2DEG 
(MODFET) wafer with the mobility and electron den- 
sity ^~70 m 2 /Vs and ?i c =2.0xl0 15 m~ 2 at 4.2 K, re- 
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spectively. The structure of the wafer was (from the 
front surface) 10 nm GaAs cap layer, 40 nm Si-doped 
(7V S i=2xl0 24 m~ 3 ) Alo.33Gao.e7As layer, 40 nm undoped 
Alo.33Gao.67As spacer layer, and 1 /xm GaAs layer with 
2DEG channel residing near the interface to the upper 
layer. The electron density, hence the mobility, was var- 
ied by illumination with infrared light-emitting diode 
(LED) when necessary. The potential modulation was 
introduced by placing a grating of high-resolution nega- 
tive electron-beam (EB) resist (calixarene) 15 on the sur- 
face of a 2DEG wafer patterned as a Hall bar by wet 
etching. The grating was oriented perpendicular to the 
direction of current so that the current flows across the 
modulation. The direction of the current was chosen to 
be a (110) direction. On cooling down the sample to 
cryogenic temperatures for measurement, elastic strain 
results from differential contraction between the resist 
and GaAs, which couples to 2DEG by deformation po- 
tential 10 or piezoelectric effect 12 to introduce modula- 
tion, the latter mechanism being maximized and dom- 
inant for our choice of the crystallographic direction. 13 
The depth d~99 nm of the 2DEG from the surface was 
evaluated as the sum of the distance between the surface 
and heterointerface, 90 nm, and the average distance of 
the 2DEG wave function from the interface, the latter 
being assessed as 8.7 nm (at n o =2.0xl0 15 m -2 ) by a 
numerical self-consistent calculation including exchange 
and correlation effect. 16 ' 17 Note that the values of a are 
kept smaller than twice the d in all of our ULSL sam- 
ples. Magnetoresistance measurements were carried out 
at either 1.4 or 4.2 K (separate cryostats were used for 
each temperature), employing a standard low- frequency 
ac lock-in technique. We assume throughout the paper 
that the state of elastic strain hence the modulation am- 
plitude do not differ between the two temperatures. 

3. Results 

Figure 1(a) shows typical magnetoresistance traces 
at 1.4 and 4.2 K, exhibiting positive magnetoresistance 
(PMR) 18 emanating from B=0 and the commensurabil- 
ity oscillation (CO), the latter being the subject of main 
interest in the present paper. In a previous paper, 7 the 
present authors have shown that CO for sinusoidal po- 
tential modulation V{x)—V\ cos(2irx/a) is well described 
by, 



osc 
J xx 



Po 



A 



|B|-^sin 



2tt- 



2R C 



where R c —hk-p/e\B\ represents the cyclotron radius with 
fcp=\/27m e the Fermi wave number, and 7 a constant 
determined by sample parameters /i and n c as 



2V2r7 \ e 



2m* 



The factors A(tt/ fiwB) and A(T/T a ,i) specified using 
the function A{x)=xj sinh(x) express damping due to 
scattering and temperature, respectively, of the oscilla- 
tion amplitude, which would otherwise decrease simply 
in proportion to the decreasing magnetic field. The pa- 
rameter is a measure of the scattering that diverts 
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Fig. 1. (a) Magnetoresistance (left axis) at 1.4 K (solid trace) and 
4.2 K (dotted trace) for ULSL with a=184 nm and n o =2.3xl0 15 
m~ 2 . The horizontal marker indicates the window to be used 
for the Fourier analysis. The damping factor A(n///,yfB) (right 
axis) due to scattering is plotted by dash-dotted line for /ivv=H 
m 2 /Vs obtained by fitting to the 1.4 K trace, (b) Oscillatory part 
of magnetoresistance obtained by subtracting a slowly-varying 
background. Calculated traces using eq. (1) are also plotted by 
thin solid (1.4 K) and dotted (4.2 K) lines, (c) Thermal damping 
factors A(T/T ,„) at 1.4 K (solid lines) and 4.2 K (dotted lines) 
for n=l, 2, and 3. See eqs. (4) and (5). 25 



electrons away from the cyclotron orbit and have been 
shown 7 to be essentially identical with the quantum mo- 
bility hq deduced from the analysis of the damping of the 
Shubnikov-de Haas (SdH) oscillation. 19 The characteris- 
tic temperature T tti i for the thermal damping factor is 
given by k^T a} i={\ /27r 2 )(afc F /2)fiw c with lu c = eB/m* 
the cyclotron angular frequency. The temperature differs 
from its counterpart in the SdH oscillation T c only by 
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the factor akp/2, and therefore the sensitivity of CO to 
temperature approaches that of SdH for small a. The de- 
pendence of A(tt I B) and A(T/T ai i) on B is plotted 
in Figs. 1(a) and 1(c), respectively. The figures clearly 
demonstrate that the damping factors are indispensable 
for the understanding the envelope of CO; particularly, 
the thermal damping factor A(T/T a ^i), 6,20 although of- 
ten neglected in attempts to explain CO amplitudes, 2 1-23 
takes on greater importance as the period a becomes 
smaller. 

As illustrated in Fig. 1(b), the oscillatory part of 
the magnetoresistance, obtained by subtracting a slowly- 
varying background following the prescription delineated 
in ref. 7, shows excellent agreement with the calculated 
trace given by eq. (1), where and V\ are used as 
fitting parameters. The agreement is especially good at 
4.2 K, where the experimental and calculated traces are 
indistinguishable in Fig. 1(b). At 1.4 K, deviation can 
be observed for B above ^0.13 T, mainly attributable 
to the disturbance by the SdH effect. It is tempting to 
take the good agreement as evidence for the modulation 
being actually simply sinusoidal. We will show below, 
however, that the agreement does not necessarily elimi- 
nate the possibility of the presence of higher harmonics 
in the modulation profile. 

Equation (1) can readily be extended for potential 
modulation including higher harmonics, 



V(x) = ^ V n cos I — nx J , 



(3) 



since each harmonics contribute independently. 24 The 
expression of CO then becomes 



,osc 

XX 



7T \ 



T \ nVj 



lA " \\B\YA — K^sin 2n 
\vwBj 1 \T a , n J a V a/n 

(4) 

The characteristic temperature for the n-th harmonic re- 
duces to 1/n of that of the fundamental component as 

1 (a/n)k F huj c T a 

where t =(1 /2-K 2 )(eh/2m*){ak? / 'fee) • The reduction in 
the characteristic temperature deeply influences the way 
the harmonic reveals itself in CO. The factor A{T/T a ^ n ) 
decreases with decreasing B from unity, the asymptotic 
value at the high-B limit, to zero; the decrease is more 
rapid for higher T, or for smaller T a /n which reflects 
either the decrease in a or the increase in n. As illus- 
trated in Fig. 1(c), the degree of damping is quite sensi- 
tive to the change in T and n in the magnetic-field and 
the temperature range of the present interest. As a re- 
sult, higher harmonics are much more heavily damped 
than the fundamental component; A(T/T a , n ) works as a 
"low-pass filter" . This makes it more difficult for the CO 
with larger n to be observed in the magnetoresistance 
traces, especially at the higher temperature. 

The presence of higher harmonics, despite the seeming 
absence in Fig. 1(b), can be testified by Fourier spectra 
shown in Fig. 2(b). For the Fourier transform to be effec- 
tive, it is advantageous to extract in advance oscillatory 
part from the raw magnetoresistance data. In Fig. 2, we 
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Fig. 2. (a) The second derivatives (d 2 /dB 2 )(Ap xx /po) of mag- 
netoresistance traces shown in Fig. 1(a), plotted against inverse 
magnetic field (within the window show by the horizontal marker 
in Fig. 1(a)). (b) Fourier spectra of the traces in (a). The spec- 
trum for 4.2 K is multiplied by a factor 3.79 to align the n=l 
peak height. 



made use of second derivative with respect to B as a 
convenient way to separate out the oscillatory part. The 
resulting (d 2 /di? 2 )(Ap a;a ;/po) is plotted against inverse 
field in Fig. 2(a) for both 1.4 and 4.2 K, whose Fourier 
transforms are the spectra shown in Fig. 2(b). The spec- 
tra exhibit peaks from the fundamental and the second 
and third harmonic CO as well as a peak from SdH. The 
spectra are normalized so as to equalize the height of the 
n=l peak for both temperatures; it can readily be ob- 
served that the peak height for T=4.2 K relative to that 
of T=IA K becomes progressively smaller for n—2 and 
3 CO, reflecting the enhanced temperature dependence 
of A(T/T a . n ) for increased n (smaller T a . n ). The SdH 
peak is not resolved for T=4.2 K for the magnetic-field 
window used in the present analysis, because of its still 
higher sensitivity to the temperature. 

It should be cautioned that the second derivative 
works sort of high-pass filter that emphasizes 

rapidly varying component, the higher harmonics in 
the present context. This remains innocuous for qual- 
itative argument of, say, the temperature dependence 
of a particular peak. For quantitative comparison be- 
tween different peaks, however, this approach is ob- 
viously inadequate. Instead, raw magnetoresistance af- 
ter subtracting the slowly-varying background as was 
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Fig. 3. Oscillatory parts of the magnetoresistance (replot of Fig. 
1(b)), plotted against inverse magnetic field (within the window 
show by the horizontal marker in Fig. 1(a)). Calculated curves 
using eq. (1) are also plotted by thin solid (1.4 K) and dotted 
(4.2 K) lines, (b) Fourier spectrum of the 1.4 K trace in (a) after 
divided by the B-dependent prefactor Fi(B). 



done in Fig. 1(b) should be analyzed directly. Accord- 
ing to eq. (4), Ap° s £/po includes the n-th harmonic 
each weighted by the B- and n-dependent envelope hmc- 
tion F n (B)=A(ir / p w B)A{T /T a , n )\B\. Therefore, a value 
proportional to nV^ is obtained by performing Fourier 
transform to Ap° s £ / po divided by F n (B), plotted as a 
function of B . An important point to keep in mind is 
the necessity to use different F n {B) depending on n to 
deduce nV^ from the same Ap° s £ / po- Figure 3(b) shows 
a Fourier spectrum taken from Ap° s £/po, shown in Fig. 
3(a), after divided by F^B). The peak at 0.82 T is pro- 
portional to 1 x Vi- Evaluation of 2V 2 2 and 3V 3 2 , how- 
ever, requires Fourier spectra different from that shown 
in Fig. 3(b), namely, the ones obtained from the oscilla- 
tory part divided by i^(-B) and F^(B), respectively, in- 
stead of Fi(B). Repeating the Fourier transforms, we can 
get the values proportional to 2V 2 2 and 3V^ 2 , hence the 
ratio IV2/V1I and IV3/V1I, from the heights of the peaks 
at 1.64 and 2.46 T of the appropriate Fourier spectra. 
Note that only the absolute value of V n can be deduced 
from CO, since the amplitude of the CO depends only on 
the square of V n 's [see eqs. (1) and (4)] and their signs 
are irrelevant. 

There is another, more practical way to evaluate n=2 
and 3 harmonic components from a single Fourier spec- 



trum such as shown in Fig. 3(b); the function to be 
Fourier transformed can be rewritten as that for n=l 
multiplied by F 1 {B)/F n {B)=A{T/T a>1 )/ A(T/T a>n ), and 
since the correction factor is much more slowly varying 
function of b = B^ 1 than the rest of the function in- 
cluding CO, it can be taken out of the integral of the 
Fourier transform and replaced by its average over the 
window [BJ ,B^~ ]=[bf,bi] for the transformation. The 
corrected peak heights are thus obtained by multiplying 

r b * A(ThlrA 

-db (6) 

'bf 



/ A(T/T aA ) 
\A(T/T a , n ) 



f 

Jbr 



b f J bf A(nTb/ Ta ) 

to the corresponding peaks in Fig. 3(b). Equation (6) can 
readily be calculated analytically for n—2 and 3. 26 

In order to determine harmonic contents \V n \ (n=l, 
2, 3) as accurately as possible from CO, we take the 
following strategy. First, we obtain \ V\\ by fitting eq. (1) 
to Ap° s £/p at 4.2 K. The higher temperature is selected 
because higher harmonics are more effectively damped 
(see Fig. 1(c)) and therefore their contribution to CO can 
be more safely neglected. Second, the ratio | V2/V1 1 and 
I V3/V1 1 are deduced from Fourier spectra using Ap° s £ / po 
at 1.4 K, since sensitivity to higher harmonics increases 
at lower temperature. The analysis is done for Ap° s £ / po 
normalized by adequate envelope function for n, F n (B), 
or employing F±(B) and then using the correction factor 
eq. (6), both giving consistent results. Finally, | V2 1 and 
I V3 1 are calculated from the value of \V\\ and the ratios 
I V2/V1 1 and I V3/V1 1 . Here we rely on our assumption that 
the harmonic contents do not vary with temperature, so 
long as the temperatures are already low enough. 

The ratios and \V n \'s thus obtained for five ULSL 
samples are plotted in Figs. 4(a) and 4(b), respectively, 
against the period a. The second harmonic is observed for 
ULSL with a>138 nm, and the third harmonic only for 
a=184 nm. Note however that the absence of the corre- 
sponding peaks in the Fourier spectra for smaller period 
ULSL does not necessarily attest to the nonexistence of 
the higher harmonics. It simply reveals that they are less 
than the limit of their observation through CO; the limit 
becomes severer for smaller a because of the dependence 
of the thermal damping factor on a, as has already been 
discussed thoroughly above. 

The ratios |U 2 /Ul| and |U 3 /Ui| plotted in Fig. 4(a) were 
acquired from the average of typically ten data points of 
the analyses described above, repeated for values of n c 
ranging from 1.9 to 3.1 xlO 15 m~ 2 varied by LED illu- 
mination. We expected, and actually observed, that the 
ratios do not show systematic dependence on n c . The 
averaging was done to reduce statistical scattering that 
could have crept in during the measurement and/or the 
analysis, and the error bars represent the standard de- 
viation. The ratios are roughly 0.3 and do not depend 
much on a. Figure 4(b) illustrates that |Vi| rapidly de- 
crease with decreasing a. The plot is for n e =2.0xl0 15 
m~ 2 , and |Vi| showed slight decrease with increasing n e 
as has been detailed in ref. 17. The dependence of \V n \ on 
a will be discussed in the subsequent section. It is easy to 
show that the oscillatory part reconstructed using eq. (4), 
substituting the values of |V^|'s obtained here, introduces 
only negligibly small modification to that calculated by 
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Fig. 4. (a) The ratio of the amplitude of the 2nd and 3rd har- 
monic potential to that of the fundamental component, obtained 
by the analysis of the Fourier spectra detailed in the text, plotted 
against the period of ULSL. (b) The amplitude of the potential 
modulation versus the period. The fundamental component Vi 
was obtained by fitting the experimental 4.2 K traces to eq. (1). 
Harmonics V2 and V3 were deduced using V\ and the ratio plot- 
ted in (a). Vi vs a is fitted to eq. (12) with ui=4.8 nm, and using 
<j>i and d as fitting parameters; excellent fitting was attained 
with 1 0i I =72.3 mV and d=97.7 nm. Dotted line for V2 is just an 
eye-guide. 



cq. (1) using |Vi| alone, owing to the V 2 dependence as 
well as to the factor A(T '/T a-n ), in full support of the 
properness of the prescription we have followed to de- 
duce |T^i|'s. 

4. Discussion 

4-1 The grating-induced effect localized at the surface 

We start by a simplistic model that (i) any alteration 
in the GaAs/AlGaAs crystal introduced directly by our 
grating is localized at the surface, and (ii) the charge that 
might be present between the surface and the hcteroint- 
erface is negligibly small. Then the electrostatic potential 
fi(x, z) inside the crystal is given by the solution of the 



Laplace's equation V 



)=0, with the boundary con- 



ditions imposed by the grating at the surface, 



fi(x, 0) = ^ <fin cos 



71=0 



2tt 



-nx 



(7) 



and deep in the substrate, (f>(x, z) — > (z — > 00). Here, we 
have taken the z-axis normal to the 2DEG plane pointing 
inward the crystal, with z=0 denoting the surface. The 



harmonic function that satisfies the boundary conditions 
is simply given by 



fi(x, z) = ^2 4>n cxp I 



n=0 



2tt 



2tt 



-nx 



(8) 



with fin's defined in cq. (7). The potential modulation 
seen by the electrons at the 2DEG plane at z=d, eq. 
(3), 27 is given after taking into account the screening 
by the 2DEG electrons themselves. Using Thomas-Fermi 
approximation for the screening, the Fourier components 
become 

V n = (-e)0„exp [ -— nd) / e TF (—n) , (9) 



with 



£tf (<?) 



(10) 



where a^=47reoeh 2 / m* e 2 represents the effective Bohr 
radius (10.4 nm for GaAs with the relative dielectric 
constant e— 13.18 and the effective mass m*=0.067m c ) 
and F(q) the form factor reflecting the finite thick- 
ness of the electron wave function. 17, 28-30 The Thomas- 
Fermi approximation is appropriate for the length scale 
longer than half the Fermi wavelength, 17, 29,31 which is 
actually the case for the range of n c encompassed in 
the present study, even for our shortest length scale 
184/3 nm. Note that exF(g) decreases with increas- 
ing a, 32 so that smaller length scale (higher harmon- 
ics) survives the screening better. The form factor may 
be evaluated numerically using numerically calculated 
wave function. 17 Instead, we employ in the present 
study F(q) calculated using Fang-Howard wave func- 
tion 33 ( FH (z)=e{Z)yfb T /2Zexv{-bZ/2) with Z=z-z 
and z representing the location of the hcterointerface 
(20 =90 nm for the present samples) and &(Z) the unit 
step function, 



l + (9/8)(q/b) + (3/8)(q/b) 2 
(l + g/6) 3 



(11) 



and e TF ™(q,b) = l + 2F FH (q)/a* B q. The analytical for- 
mula is of benefit to the fitting to the experimental plots. 
The parameter b is a measure of the thickness of the 
2DEG wave function: w=^/3/b with w the rms thick- 
ness. Instead of using the well-known variational formula 
that determines b from the depletion charge and n , 31 
we adopt the value of 6 that equalizes w to the value 
deduced by numerically solving the Shrodinger equation 
and Poisson's equation self-consistently, as described in 
ref. 17. The formula for the fundamental component, 



Vi = {-e)fii exp 




^fe^ , (12) 
\ a w I 



with u>=4.8 nm (the value for n c =2.0x 10 15 m~ 2 ) is fitted 
to the experimental plot using fix and d as fitting param- 
eters, and is also shown in Fig. 4(b). The excellent agree- 
ment is achieved with the parameters |</>i|=72.3 mV and 
d=97.7 nm, the latter being quite close to the designed 
value of 98.7 nm. The value of fi\ is also reasonable, as 
will be discussed later. The good agreement suggests that 



6 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



the model adopted to deduce eq. (12) gives a satisfactory 
description for the fundamental component, namely, the 
origin of V\ can be traced back to be located at the sur- 
face. For V2 and V3, similar fitting is difficult because of 
the lack of enough number of data points. However, even 
a cursory glance at the plot reveals it difficult to describe 
V2 well with eq. (9): according to eq. (9), V2 should de- 
crease more rapidly with decreasing a than V\, which is 
the opposite to the situation observed in Fig. 4(b). So 
far, we have tacitly assumed that the </>„'s determined 
by the boundary condition at the surface do not depend 
on a, a reasonable assumption since we prepared differ- 
ent ULSL samples under the same conditions, employing 
the same GaAs/AlGaAs 2DEG wafer and the identical 
combination of EB resist and dose, hence having prac- 
tically the same resist-film thickness (~30 nm) and the 
elastic properties. The dependence of V n on a will alter 
if we discard this assumption. However, it is still quite 
unlikely that the observed V2 and V3 are compatible with 
eq. (9). Their ratio to V\ using eq. (9), 



Vn 
V 1 



0„e T F(27r/a) 



exp 



- — (n-l)d 
a 



(13) 



6i£TF(27rn/a) 

becomes very small because of the exponential factor, un- 
less the ratio 4> n /4>i is unacceptably large. For example, 
\4>2/4>\\— 5 and |</>3/</>i|~100 is required to explain the 
experimental | V2/V1 j and | V3/V1 1 by eq. (13) for o=184 
nm, and the values are still larger for smaller a. This 
implies that the origin of the second and the third har- 
monics must be sought after in the grating-induced effect 
beyond the effect localized at the surface. 

4-2 Comparison with the theories of strain-induced po- 
tential modulation 
The strain induced by the grating is, of course, not 
perfectly localized at the surface but propagate in- 
side the crystal. Therefore the strain can in principle 
generate potential modulation at an arbitrary depth 
z through deformation potential or piezoelectric ef- 
fect. Analytic expressions are presented in the theo- 
ries of strain-induced potential modulation. 10, 12 Here 
we quickly review the results of the theories essential 
to the present discussion. Writing Airy's stress func- 
tion as z\ with x a harmonic function, xi x > z ) = 
^7=0 Xn exp (-q n z) cos (q n x), with Xn's determined by 
elastic boundary condition (compatibility with the grat- 
ing) at the surface x(x, 0) and q n =2nn/a, Fourier coef- 
ficients for the deformation potential is given by, 

S 2(1 + i/)(1-2i/) 



i-e) 



E 



-QnXn exp (-q n z) 



(14) 

where H, v, and E represent the deformation potential 
constant, Poisson's ratio, and Young's modulus of the 
host semiconductor crystal, respectively. The piezoelec- 
tric charge density p(x, z) for the angle 0=45° between 
the [100] axis and the x axis reads in Fourier coefficients 

Pn{z) = ^14 [(5 - 2u) - 3q n z] q n Xn exp (-q n z) , (15) 

with du the only surviving element of the piezoelectric 
tensor in the crystal with F43m symmetry. The piezo- 



electric potential is found by solving Poisson's equation 
V 2 tp pz (x, z)~—p(x, z)/e e. It can generally be written as 

Vn\A = [Cn + (7 - A V )q n Z - 3(qn z ) 2 ] Xn CX P (-Q n z) , 

(16) 

with C n constants to be determined by boundary condi- 
tions. The condition ?/> pz (x, z) — > (z — > 00) is automat- 
ically satisfied. In rcf. 12, an cquipotcntial layer at z~c 
(either the surface or the 5-doped layer), ^> pz (a;,c)=0, j s 
considered, which results in C„=3(g„c) 2 — (7 — Av)q n c. 
For electrons at 2DEG, V n icf =(-e)ij n cf (d)/e TF (q n ) and 
VP z =(-e)ii n z (d)/e TF (q n ) after screened by the 2DEG. 35 
Note that both ip n cl (z) and ip%?(z), as well as 
(f> n exp>(—q n z) in the previous section, include the factor 
exp(— q n z) and therefore attenuate with z. The attenua- 
tion is more rapid for smaller a or for higher harmonics. 
What is peculiar in ip pz (z) is its inclusion of the poly- 
nomial of q n z, leading to non-monotonic behavior with 
z around z^a/n. Therefore, V pz /V pz can become large 
when z—d happens to lie near the node of ip pz (z) without 
resorting to unnaturally large Xn/xi- 

Unfortunately, it is not possible to make a quantitative 
comparison between our data and the theories, since the 
evaluation of x„'s requires the knowledge of the elastic 
constants for our EB resist, which are currently unavail- 
able. 36 For qualitative comparison of the dependence on 
a, we plot V^ ci /xn and V pz /x n (n=l, 2, 3) against a in 
Fig. 5, with eTF(9n)=eTF FH (<Zn, y/3/w) and using elas- 
tic parameters for GaAs 12,37 and structure parameters 
for our present samples. We also plot, with a separate 
(right) vertical axis, V\/{— e<f>i) calculated by eq. (12), 
which describes the experimental behavior well, and its 
equivalent for n=2 and 3. For the calculation of V£ z , we 
adopted two different models for boundary conditions at 
the surface after refs. 10 and 12. In the "pinned" sur- 
face model, appropriate if the sample is in full equilib- 
rium, we assume the surface z=c=0 to be the equipo- 
tential plane tp% z (x, c)=0. 38 The "frozen" surface model 
assumes that charges at the surface state are frozen be- 
cause of slow equilibration rate at low temperatures be- 
tween the surface state and the 2DEG. Capacitance spec- 
troscopy 39 and n c vs back-gate voltage measurements 40 
suggest that the frozen surface model is actually a better 
description at cryogenic temperatures. Appealing to the 
high dielectric constant of the semiconductor, we may 
take dip n z (x,z)/dz=0 at z=0 in this model, leading to 
C„=7 — Av in eq. (16). In contrast to the situation consid- 
ered in ref. 10, we take dip n z (x, z) /dz~0 not only at open 
surface between the resist gates but also under the gates, 
because our grating is made of non-conducting material; 
we do not need to solve mixed boundary value problem. 
For both boundary conditions, Vf z shows a maximum 
in the range of a investigated in the present study, at 
variance with the experimental result. The dependence 
of Vf cl on a appears to be closer to that of the exper- 
iment. However, as mentioned earlier and is well estab- 
lished, 13 the strain-induced effect is dominated by the 
piezoelectric effect, V pz ^>V^ ei , for the present crystallo- 
graphic direction, which is actually observed in Fig. 5. 
Therefore, neither V 1 pz nor gives better account of 



jll Paper 



Author Name 7 



0.30 
0.25 
0.20 
0.15 
0.10 
0.05 
0.00 



rt=1 




piezoelectric . -'* 
pinned 




frozen 




deformation, 

// x10,,>' 






— surface 








80 100 120 140 160 180 200 
a (nm) 



Fig. 5. Normalized piezoelectric and deformation potentials (left 
axis), V£ z /xn and /Xn (n=l, 2, 3), as a function 

of a, calculated with eqs. (16) and (14), respectively, and 
eTF(<7«)=eTF FH (<7n, VS/w), for our sample parameters d=98.7 
nm, w=4.8 nm. Elastic parameters for GaAs, 12 ' 37 S=— 8.2 eV, 
^=0.31, £=85.3 GPa, and di 4 =-2.69x lCT 12 mV" 1 arc used. 
For V„ z /xn, two different boundary conditions at the surface, 
"pinned" and "frozen" surface models, are considered and are 
plotted by dotted and solid lines, respectively. V„ /xn (dashed 
lines) arc magnified by a factor 10. V n /(—erf> n ) calculated by eq. 
(9) with £tf (Qn)=£TF FIi (in , V3/w) are plotted by dash-dotted 
line (right axis) for comparison. Note the expanded (x5) vertical 
scale (both right and left axes) for n=3. 



the experimental V\, compared with the success of the 
model described in §4.1. An attempt to explain this will 
be presented in the next section. For n—2 and 3, our 
experimental data arc neither sufficient in quantity nor 
accurate enough for detailed comparison with the the- 
ory to be meaningful. Nonetheless, Fig. 5 does suggest 
that the piezoelectric effect is the most likely origin of 
the experimentally observed large |V^/Vi| {n—2, 3). 

4-3 Plausible model for our potential modulation 

In the theories, 10, 12 several simplifications are made to 
enable analytical treatments. Above all, elastic proper- 
ties of the layers between the surface and the heteroint- 
erface are all approximated by those of a single material, 
GaAs, and are treated as isotropic. We suspect this sim- 
plification to be partly responsible for the disagreement 
between the behavior of our data and eq. (16). In re- 
ality, 80 nm out of the 90 nm between the surface and 
the heterointerface is AlGaAs, with the upper half of the 
AlGaAs layer doped with Si. Slightly different clastic, 



thermal, or electrostatic properties of each layer stacked 
successively can in principle alter the potential intro- 
duced to the 2DEG through the layers. Specifically, the 
lattice mismatch between GaAs and AlGaAs is about 
0.1%, roughly the same order of magnitude as the esti- 
mated strain introduced by the grating via differential 
contraction. 10, 12 Therefore the stress introduced to the 
GaAs cap layer at the surface from the resist is possibly 
modified when transmitting across the interface, 10 nm 
below the surface, to the Si-doped AlGaAs layer; depend- 
ing on the sign of the strain at the surface, the stress will 
be either enhanced or reduced at the interface. If the sign 
and the magnitude of the strain is such that the stress is 
nullified at the interface, the piezoelectric charge and po- 
tential may taken to be generated only in the vicinity of 
the surface, realizing the situation assumed in §4.1. We 
speculate this is basically what is taking place for the 
fundamental (n=l) component. We take here a further 
step to see whether the value of <j>i obtained from the 
experimental data in §4.1 is reasonably explained within 
this model by roughly estimating the size of the differ- 
ential contraction required to obtain the value. 

For simplicity we assume the stress hence the piezo- 
electric charge to be generated only at the surface, ne- 
glecting the finite thickness of the cap layer along whose 
depth the stress is expected to attenuate. We thus use eq. 
(16) with z=0 for the piezoelectric potential that serves 
as a boundary condition eq. (7) for the electrostatic prob- 
lem, resulting in 



— C y 

8e e 



(17) 



The relation of \{ x i 0) with the differential contraction 
a gate AT (agate denoting the difference of the linear ex- 
pansion coefficient between the gate and the semiconduc- 
tor) is given by the compatibility equation (the bound- 
ary condition for the elastic problem at the semiconduc- 
tor surface). In an approximation named "elastic" gate 
model in rcf. 12, the compatibility equation reads 



X ( X ,0)= hE T ea ^ eAT u(x), 



1 



^gatc 



(18) 



where h, t'gate, and i? ga te represent the thickness, Pois- 
son's ratio, and Young's modulus for the gate, and 
u(x)=l for x under the gates and for x in between. The 
"elastic" gate model applies to thin gates. There the gate 
is forced to contract with the same ratio as the under- 
lying semiconductor, uniformly stressed; the gate reacts 
back a force, concentrated at the edge, on the semicon- 
ductor. The model appears to be suitable for our thin 
(fr~30 nm), and probably soft, resist gate. For a grat- 
ing with duty ratio a (^1/2 in the present samples), 
u{x) = a+ (2/ir) X)^Li[ sm ( a7rn )/ n ] cos (Qnx), thus eq. 
(18) becomes in Fourier coefficient 

/legato a g ato AT sin(a7m) 



Xr, 



1 



(19) 



^gate 7Tn / 2 

for n=l, 2, 3,.... From eq. (17) and using the "frozen" 
surface model, C„=7— 4^, \xi | =4.4 N/m is obtained from 
|0i |=72.3 mV. Equation (19) shows that the value of \i 
can be achieved with a ga t e AT considerably smaller than 
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1% by using reasonable values for elastic constants for the 
resist (e.g., ^gatc^O.3 and E gatc tens of GPa), although 
their exact values are unavailable. 

Of course it is highly unlikely that the stress intro- 
duced at the surface completely relaxes at the interface 
between the cap and the Si-doped layer. The stress and 
therefore the piezoelectric charge probably remain be- 
low the interface, although smaller than predicted by eq. 
(15). The discussion we have presented so far suggests 
that for the fundamental (n=l) component the resid- 
ual piezoelectric charge can be taken to be negligibly 
small. For higher harmonics, the remnant piezoelectric 
charge takes on more important role, since the strong 
piezoelectric potential limited to the vicinity of the sur- 
face becomes less influential at the 2DEG plane due to 
the exponential factor, as discussed in §4.1. We presume 
that the residual piezoelectric charge is mainly respon- 
sible for the higher-harmonic potential modulation seen 
by the 2DEG electrons. An alternative candidate for the 
origin of the higher harmonics are the charge building 
up at the heterointerface just above the 2DEG due to 
the difference in the elastic, piezoelectric, or dielectric 
constants, as discussed in ref. 12. It seems quite difficult, 
at this stage, to further elaborate and specify the origin 
from our currently available data. 

5. Conclusions 

We have investigated the harmonic content included 
in strain-induced potential modulation, with recourse to 
detailed analysis of CO in ULSL's. We have pointed out 
that the magnitude of higher harmonics can be larger 
than seemingly suggested by their appearance in raw 
CO traces, since the CO is heavily weighted in favor 
of slower oscillation due to the thermal damping factor 
A(T/T aiTl ). In our present samples, | V2 1 and | V3 1 are as 
large as ~0.3|Vi|, although the second or third harmonics 
are not readily discernible in the raw magnetoresistance 
data. 

Higher harmonics V n (n>l) have a general tendency to 
be suppressed relative to the fundamental component V\ 
owing to the ubiquitous exponential factor, exp(— q n z), 
common to all the mechanisms discussed in §4. How- 
ever, for the piezoelectric effect, the dominant effect in 
the present samples, the dependence of V n on q n z be- 
comes complicated because of the extra polynomial de- 
pendence, which we take as one of the causes of rather 
large ratio V n /V\ (n=2,3). The dependence of V n on a 
is naturally explained by assuming that the source of V\ 
is localized at the surface with negligible charge between 
the surface and the 2DEG, while V2 and V3 originate 
from piezoelectric (or interface) charge located closer to 
2DEG. Noteworthily, the screening by 2DEG electrons 
favors higher harmonics, giving advantage of the factor 
of roughly 2 and 3 for V2 and V3 , respectively. 

We want to stress that the higher harmonics are shown 
to be able to introduce into the 2DEG plane the length 
scale substantially smaller than its depth d~99 nm from 
the surface, the smallest of those detected by CO in the 
present experiments being a/2=69 nm and a/3=61 nm, 
respectively. This suggests a possibility of realizing ULSL 
with the period smaller than d by exploiting the strain- 



induced piezoelectric effect, using metallic grating; since 
the modulation introduced electrostatically by applying 
a bias to the metallic grating is almost exclusively com- 
posed of the fundamental component owing to the factor 
cxp(— q n z), it is possible to tune the bias to cancel out 
the fundamental component due to the strain-induced 
effect, thereby making the higher harmonics selectively 
survive. 41 Lateral superlattice samples with small period, 
close to or shorter than, e.g., the Fermi wave length, 
made from a minimally disordered 2DEG residing at 
large depth in the MODFET structure, will prove in- 
valuable to experimental studies of many aspects of su- 
perlattices. 
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